Astronomy & Astrophysics manuscript no. bgemoscpaper 


© ESO 2013 


March 4, 2013 





The Mass of the Planet-hosting Giant Star f3 Geminorum 
Determined from its p-mode Oscillation Spectrum 

A. P. Hatzes 1 , M. Zechmeister 2 , J. Matthews 3 , R. Kuschnig 4 , 
G.A.H. Walker 5 , M. Dollinger 1 , D.B. Guenther 6 , A.F.J. Moffat 7 , S.M. Rucinski 8 , and D. Sasselov 9 

1 Thuringer Landessternwarte Tautenburg, Sternwarte 5, D-07778 Tautenburg, Germany 

2 Georg-August-Universitat Gottingen, Friedrich-Hund-Platz 1, 37077 

3 Department of Physics, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada 

4 Institut fiir Astronomie, Universitat Wien Tiirkenschanzstrasse 17, A-1180 Vienna, Austria 

5 1234 Hewlett Place, Victoria, BC V8S 4P7, Canada 

6 Institute for Computational Astrophysics, Department of Astronomy and Physics, Saint Mary's University, Halifax, 
NS B3H 3C3, Canada 

7 Observatoire Astronomque du Mont Megantic, Department de Physique, Universite de Montreal C. P. 6128, 
Succursale : Centre- Ville, Montreal, QC H3C 3J7, Canada 

8 Department of Astronomy and Astrophysics, University of Toronto, Toronto, ON M5S 3H4, Canada 

9 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 102138 

Received; accepted 

ABSTRACT 

Aims. Our aim is to use precise radial velocity measurements and photometric data to derive the frequency spacing of 
the p-mode oscillation spectrum of the planet-hosting star ft Gem. This spacing along with the interferometric radius 
for this star can then be used to derive an accurate stellar mass. 

Methods. We use a long time series of over 60 hours of precise stellar radial velocity measurements of ft Gem taken 
with an iodine absorption cell at the echelle spectrograph mounted on the 2m Alfred Jensch Telescope. We also present 
complementary photometric data for this star taken with the MOST microsatellite spanning 3.6 d. A Fourier analysis 
is used to derive the frequencies that are present in each data set. 

Results. The Fourier analysis of the radial velocity data reveals the presence of up to 17 significant pulsation modes in 
the frequency interval 10-250 ^iHz. Most of these fall on a grid of equally-spaced frequencies having a separation of 7.14 
± 0.12 /iHz. An analysis of 3.6 days of high precision photometry taken with the MOST space telescopes shows the 
presence of up to 16 modes, six of which are consistent with modes found in the spectral (radial velocity) data. This 
frequency spacing is consistent with high overtone radial pulsations; however, until the pulsation modes are identified 
we cannot be sure if some of these are nonradial modes or even mixed modes. The radial velocity frequency spacing 
along with angular diameter measurements of ft Gem via interferometry results in a stellar mass of M = 1.91 ± 0.09 
Mq. This value confirms the intermediate mass of the star determined using stellar evolutionary tracks. 
Conclusions, ft Gem is confirmed to be an intermediate mass star. Stellar pulsations in giant stars along with inter- 
ferometric radius measurements can provide accurate determinations of the stellar mass of planet hosting giant stars. 
These can also be used to calibrate stellar evolutionary tracks. 

Key words, star: individual: ft Gem, - techniques: radial velocities - stars: late-type - stellar oscillations 



1. Introduction 

There are several programs underway that use precise stel- 
lar radial velocity (RV) measurements to search for planets 
around evolved giant stars (e.g. Setiawan et al. 2003; Sato et 
al. 2003, Johnson et al. 2007; Niedzielski et al. 2009; Han et 
al. 2010). The reason is that the progenitors of these giant 
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stars are A-F dwarfs - stars that have high effective temper- 
atures and rapid rotation rates that makes it almost impos- 
sible to detect the subtle reflex motion from any planetary 
companion with RV measurements. Although giant stars 
offer us a means of detecting planetary companions around 
stars more massive than the sun, there is a drawback in 
the determination of an accurate stellar mass. Unlike for 
main sequence stars that have a tight and reasonably well 
calibrated relationship between stellar mass and effective 
temperature, the evolutionary tracks of giants stars cover- 
ing a broad range of stellar masses all converge to the same 
region of the H-R diagram. One must rely on evolution- 
ary tracks, which are model dependent and require accu- 
rate stellar parameters (luminosity, effective temperature, 
abundance, etc.), to determine the stellar mass. Accurate 
stellar masses of planet hosting giant stars are needed if 
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Fig. 1. Radial velocity measurements for (3 Gem on 9 different nights. The error bar in each panel shows the typical error 
for each night. The solid line represents the fit to the data using the frequencies and amplitudes in Table 2. For clarity 
we used variable limits for the ordinate. 

we are to understand the process of planet formation as a 
function of stellar mass. 

The bright star (3 Gem is a planet-hosting K giant star 
(Hatzes & Cochran 1993; Hatzes et al. 2006; Refiert et al. 
2006) that was recently shown to exhibit stellar oscillations. 
Hatzes & Zcchmcistcr (2007; hereafter HZ07) using 20 hr of 
precise stellar radial velocity measurements demonstrated 
that this star shows up to 6 oscillation modes in the fre- 
quency range v = 10 - 150 /xHz. The most dominant mode 
was at v — 89 ^Hz and had an amplitude of 5 ms _1 . The 



pulsation frequencies and amplitudes were consistent with 
those expected for solar-like p-mode oscillations for a giant 
star with radius R = 8.8 Rq and mass, M w 2 M Q . The 
former was determined from interferometric measurements 
(Nordgren et al. 2001). 

The mass estimate of /? Gem by HZ07 was based on the 
frequency of maximum power, i/ mal , and the scaling rela- 
tions of Kjeldsen & Bedding (1995). This is not as accurate 
for two reasons. First, using a short time span of data (3 
nights in the case of HZ07) we cannot be sure that the dom- 
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inant mode we detect accurately represents v max . Likewise 
fitting a Gaussian to the envelope of observed excess power 
may have a large error. Second, to derive the mass from 
Vmax one requires a knowledge of the effective stellar tem- 
perature which is not well constrained. On the other hand, 
the stellar mass determined from the frequency spacing de- 
pends only on the stellar radius (measured with interfer- 
ometry in this case) and is thus independent of other more 
poorly determined stellar parameters like temperature or 
surface gravity. If one knows the stellar radius, one can thus 
derive the stellar mass. Just as important, determining the 
precise frequencies of as many detected modes as possible 
and their separation are more valuable for asteroseismic 
modeling than just determining v max Unfortunately, the 
data set of HZ07 was of insufficient length to do this. 

In this paper we use a much larger data set with over 
60 hr of coverage, including the measurements presented 
in HZ07, to derive the oscillation spectrum and thus the 
frequency spacing of (3 Gem. We also present photometric 
observations taken with the MOST space telescope. 

2. Observational Data 

2.1. Spectral Data 

Radial velocity measurements were made using the coude 
echelle spectrograph of the 2m Alfred Jensch Telescope 
of the Thiiringer Landessternwarte Tautenburg (Thuringia 
State Observatory). Observations were made on 10 nights 
over the time span 24 December 2006 to 18 April 2007, 
or a difference of 116 days from first to last observation. 
A total of over 60 hrs of observations were made result- 
ing in over 762 measurements. Precise stellar radial veloci- 
ties were determined using an iodine absorption cell placed 
in the optical path of the telescope. A further description 
of the instrumental setup and RV analysis method can be 
found in Hatzes et al. (2005). Table 1 lists the journal of 
observations. 
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Fig. 2. Residual radial velocities after subtracting the fit 
using the frequencies and amplitudes in Table [21 In each 
panel the rms scatter for that night is given. 

compact figure for easy viewing we do not show the last 
night which was also included in our analysis. Although the 
time span of the observations (time between first and last 
data point) for this night is large compared to the other 
nights, there were large gaps and only 39 measurements 
were made due to poor observing conditions. Also, because 
of these conditions the RV precision was about a factor of 
2 worse than the other nights. If we exclude these measure- 
ments the following results are not significantly altered. 
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6.52 


39 


3.13 



Table 1. The journal of observations including Julian Date 
as the start of the night, time coverage per night, number 
of observations each night, and the rms scatter of the data 
about the multi-component frequency fit (see text). 



The planetary companion to /3 Gem causes a reflex mo- 
tion of 39 ms^ 1 with a period of 590 days. This orbital 
motion was removed and the subsequent analysis performed 
on the residual RV data. Figure Q] shows the RV measure- 
ments for 9 of the nights after removal of the orbital mo- 
tion. Figure [5] shows the residual RVs after subtracting the 
fit discussed in more detail below. In order to make a more 



2.2. Photometric Data 

Photometric measurements were taken using the 
Microvariability and Oscillations of Stars Microsatellite 
(MOST) space mission as part of a secondary observational 
program. MOST is 15/17.3 cm Rumak-Maksutov telescope 
with a CCD photometer in a circular Sun-synchronous 
polar orbit with a period of 101.41 minutes (Walker et 
al. 2003). MOST observed /3 Gem on 5 consecutive nights 
spanning 13 - 17 February 2007 (JD = 2,454,145.633 and 
2,454,149.26). Individual exposures were set to 1.5 sees 
and 14 images were stacked on board before downloading 
the data. Observations of (3 Gem alternated with those of 
a primary program of MOST. /? Gem was continuously 
observed for 40 min followed by a gap of approximately 
1 hour for a duty cycle of about 40%. A total of 4440 
photometric measurements were made of /3 Gem 

Figure [3] shows the time series of the photometric mea- 
surements. For clarity data points shown in the figure are 
average values of 4 measurements (50 sees). 

3. Stellar Parameters 

(3 Gem is a nearby K0 III star at a Hipparcos-measured 
distance of 10.3 pes (van Leeuwen 2007). The most recent 
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Fig. 3. (Top) Photometry of /3 Gem taken with the MOST 
micro-satellite. The solid line represents the fit to the data 
using the frequencies and amplitudes in Table 4. (Bottom) 
Residuals after subtracting the fit. Units are in parts per 
thousand (PPT). 



interferometric measurements for this star yielded an an- 
gular diameter of 7.96 ± 0.09 mas (Nordgren et al. 2001) 
which corresponds to a radius of 8.8 ± 0.1 Rq. 

This star has also been the subject of several spectral in- 
vestigations. McWilliam (1990) measured an effective tem- 
perature, T cS , of 4850 K, a metalicity of [Fe/H] = -0.07, 
and a surface gravity, log g = 2.96. Gray et al. (2003) de- 
rived the same effective temperature, but a lower surface 
gravity (log g — 2.52), and higher metalicity ([Fe/H] = 
0.08). Allende Prieto et al. (2004) determined T cB = 4666 
± 95 K, a surface gravity of log g — 2.685 ± 0.09 and a 
metalicity [Fe/H] = 0.19. 

We also derived the stellar parameters for /? Gem us- 
ing a high signal-to-noise spectrum taken with the coude 
echelle spectrograph of the 2m telescope at Tautenburg. 
Our analysis yielded T eS = 4835 ± 50 K, log g = 2.70 ± 
0.10, and [Fe/H] = -0.07 ± 0.05. These values are consis- 
tent with the McWilliam and Gray et al. values, and are in 
good agreement with the Allende Prieto et al. values with 
the exception that our analysis does not indicate such a 
high metalicity. 

The stellar parameters of mass, radius, and age 
were det ermined using the online t ool from Girardi 
(http://stev.oapd.inaf.it/cgi-bin/param). This tool uses a 
a library of theoretical isochrones (Girardi et al. 2000) and 
a modified version of the Baysian estimation method imple- 
mented by J^rgensen & Lindegren (2005). A detailed de- 
scription of this method is given by da Silva et al. (2006). 
As input one needs the visual magnitude V, the parallax 
7T, the effective temperature T c g, and the metalicity [Fe/H]. 
For the T e g and the metalicity we used our derived values. 
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Fig. 4. (Top) The amplitude spectrum for the RV data. 
(Bottom) the window function. 
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Table 2. Pulsation frequencies from the RV data. 



The tool yields a stellar mass of M = 1.96 ± 0.19 M , an 
age of 1.19 ± 0.3 Gyr, and a radius of R = 8.29 ± 0.35 
Rq . The radius is in very good agreement with the one de- 
rived from interferometric measurements. The stellar mass 
is also consistent with the frequency of maximum power for 
the oscillations observed by HZ07. 

4. Frequency Analysis 

4.1. The RV data 

A frequency analysis was done on the full RV data set us- 
ing the program Period04 (Lenz & Breger 2004) which per- 
forms a discrete Fourier transform on the data and then 
fits the time series with the found frequencies. We are aware 
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Fig. 5. (Top) The amplitude spectrum of the 17-component 
fit to the data that was sampled in the same manner as the 
data. Random noise with a = 2 ms _1 was also added to 
the synthetic data. (Middle) The amplitude spectrum of the 
residual RVs after subtracting the fit to the data. (Bottom) 
The amplitude spectrum of random noise with a = 2 ms _1 
that was sampled in the same manner as the data. 
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Fig. 6. Phase diagrams showing the residual RV measure- 
ments phase to the frequencies found in the Fourier anal- 
ysis. For each panel the contribution of all other frequen- 
cies were subtracted prior to phasing the data. Data were 
phased in phase bins of 0.05. The gap in panel /„g is an 
artifact of the sampling and frequency. 
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Fig. 7. (Top) DFT amplitude spectrum of the MOST pho- 
tometric data. Units are in parts per thousand (PPT). 
(Bottom) the window function. The peak at 14.2 cd _1 seen 
in the window function is due to gaps in the data (see text). 

that the oscillations in a giant star may not be coherent over 
the long time span of our observations. However, as a first 
analysis we wanted to see if a multi-component since func- 
tion could adequately fit the data, and more importantly if 
the derived frequencies had a physical interpretation. The 
time unit for our RV analysis program is days, therefore 
frequency is expressed in cd _1 . In tables we also give the 
frequency in the more common units of /xHz. Figure HI shows 
the discrete Fourier transform amplitude spectra for the full 
data set as well as the spectral window. There appears to be 
excess power in two frequency ranges: 0.5 - 3 cd _1 and 6 - 
10 c d . The highest peak in the first frequency range is at 
0.86 cd _1 and for the second range at 7.54 cd _1 . The lat- 
ter corresponds to the dominant frequency found by HZ07 
based on an analysis of a subset of the data presented here. 

In order to search for multiple modes we used the stan- 
dard pre-whitening procedure. A Fourier transform of the 
data yielded the highest peak in the interval < v < 30 
cd _1 . A least squares sine fit to the data was made with 
Period04 using this frequency and its contribution was sub- 
tracted from the data. A Fourier analysis was then per- 
formed on the residuals to find the next dominant frequency 
in the data. This process was continued until no additional 
significant peaks were found. Finally, Period04 was used 
to find a solution by fitting all sine components simultane- 
ously. 

In all cases the highest peak in the amplitude spectra 
was chosen at each step of the pre-whitening procedure. We 
are aware that the spectral window, noise, finite lifetime of 
pulsation modes, etc. can cause spectral leakage and the 
identification of a spurious pulsation mode. Rather than 
making a priori judgments as to which peak is real based 
on a preliminary frequency spacing of modes (see below) we 
simply chose the strongest peak. This may result in some of 
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our identified modes being an alias frequency or the pres- 
ence of artifact frequencies due to finite mode lifetimes. 

Our analysis found a total of 17 significant frequencies 
in the RV time series. These are listed in Table 2 (denoted 
fvi - fvVr) in the order in which they were found by the 
pre- whitening process. The line in Figure Q] shows the fit to 
the RV data using these 17 modes and Figure [5] shows the 
residual RV variations after subtracting the 17-component 
fit. The final rms scatter of the data about this fit is 1.9 
ms _1 for the full data set. The rms scatter of the data 
about the final fit for each night is listed in Table 1. In their 
analysis of the nights 6-8 HZ07 obtained slightly better fits 
with rms scatter of 1.2, 1.5, and 1.7 ms" 1 , for each night 
respectively. The slightly poorer rms scatter when fitting 
the entire data set may be due to amplitude variations in 
the different modes over the time span of the measurements 
(116 days) or due to mode lifetimes shorter than the full 
time span of our observations. If the time scales of these 
possible amplitude variations are of ~ days, then fitting 
subsets of the data should result in better fits. The white 
noise level was estimated using PeriodO^ in the frequency 
range v = 10-30 /iHz (0.86 - 2.59 cd _1 ) and after subtract- 
ing the contribution of the found modes. This level was f=a 
0.15 ms" 1 . 

The statistical significance of the periods were assessed 
using a "bootstrap randomization technique" (see Kiirster 
et al. 1999). After removing the previous dominant fre- 
quency the RV values were randomly shuffled keeping the 
times fixed and a Scargle periodogram calculated (Scargle 
1982). A Scargle periodogram was used since the power is 
a measure of the statistical significance of a signal. After a 
large number of shuffles (100,000) the fraction of random 
periodograms having Scargle power greater than the data 
periodogram gave an estimate of the false alarm probability 
(FAP), that is the probability that a signal was due purely 
to noise possibly in combination with the spectral window. 
The FAPs are also listed in Table 2. FAPs with upper limits 
indicate that there was no instance in 100,000 random data 
sets where the power of the Scargle periodogram exceeded 
the data periodogram. The first 12 frequencies found by 
our pre-whitening procedure are highly significant having 
FAP < 10~ 5 . Frequencies /13 - /15 are also very signifi- 
cant having FAPs well below 1%. The final two frequencies 
are moderately significant with FAP w 2% however, since 
these fall near the derived frequency spacing (see below) 
they may indeed be real. 

We compared the amplitude spectrum of the 17- 
component fit to that of the observed amplitude spectrum. 
The top panel in Figure [5] shows the amplitude spectrum 
of the 17-component fit to data using the frequencies and 
amplitudes of Table 2 and sampled in the same manner as 
the data. Random noise with (T = 2ms _1 was also added to 
the data. The middle panel of Figure [5] shows the amplitude 
spectrum of the final residuals of the data after removing 
the 17-component fit. The bottom panel in the figure shows 
the amplitude spectrum of pure noise with a standard de- 
viation of 2.0 ms _1 that was sampled in the same manner 
as the data. 

There are several things to note about this figure. First, 
the synthetic fit with noise has an amplitude spectrum that 
is similar to the data spectrum. A 17-component sine func- 
tion with the proper sampling and noise can adequately 
reproduce most of the features seen in the amplitude spec- 
trum. Second, the amplitude spectrum of the final residuals 



HZ07 frequency 



This Work 
(MHz) 



Comment 



29.75 ± 0.54 
48.41 ± 1.06 
79.47 ± 0.76 
86.91 ± 0.37 
104.40 ± 0.49 
149.25 ± 0.63 
193.63 ± 3.07 



31.20 
36.47 
80.49 
87.24 
115.66 
144.61 
194.72 



1-day alias 



1-day alias 
1-day alias + Afo 



Table 3. Comparison of frequencies found in this work and 
HZ07. 



also looks similar to the spectrum of random noise. Finally, 
much of the structure seen in the data amplitude spectrum 
at 0.3-3 cd -1 and 6-10 cd _1 (top panel Fig. [4]) cannot be 
due to white noise combined with the spectral window. 

We now compare the results of our frequency analy- 
sis to those of HZ07 whose analysis was based only on a 
three-night subset of the data presented here. HZ07 found 
evidence for seven modes in the data and these are listed in 
Table 3. This table also lists the nearest frequencies found 
by the current analysis. Four of the modes found by an 
analysis of the full data set coincide within the errors to 
peaks found by HZ07. Two other modes (48.41 /iHz and 
36.47 /iHz) found by HZ07 are most likely 1-day aliases of 
the true frequency. The 149.3 /iHz) mode found by HZ07 
may be an alias of 137.68 /iHz. This mode differs by 6.8 
/iHz from a mode found by the current analysis. This dif- 
ference is close to the frequency spacing, Aj/, found in this 
work (see below). In conclusion, the frequencies found by 
HZ07 on the limited subset of data are for the most part 
consistent with those found by an analysis of the full data 
set. 

Figure [6] shows the phase diagram of the RV variations 
from the individual modes found in our frequency analy- 
sis. For each displayed frequency the contribution of the 
other modes was first removed before phasing the data. For 
clarity we display phase-averaged values using bins of 0.05 
in phase. Note that the phase gap in the panel for f v $ is 
an artifact of the sampling and phasing the data to this 
frequency. 



4.2. MOST Photometry 

A period analysis was also performed on the MOST pho- 
tometry. Figure [7] shows the discrete Fourier transform of 
the MOST data and the spectral window. The white noise 
in the DFT estimated at frequencies above 30 cd -1 was 
about 0.025 parts per thousand (PPT). The pre-whitening 
procedure was applied until all frequencies with amplitudes 
greater than four times the noise level were found. These 
correspond to false alarm probabilities less than 0.01. A to- 
tal of 16 frequencies were found which are listed in Table 4 
(denoted f p i - f p ie). Also listed are the false alarm proba- 
bilities calculated using the same procedure as for the RV 
measurements and 100,000 shuffles of the data. The final 
RMS scatter of the data about the fit is 0.13 PPT. 

At least six of the photometric frequencies have counter- 
parts in the RV data. The highest amplitude photometric 
modes at 2.11 and 0.38 cd -1 appear to have no counter- 
parts in the RV data. However, we caution the reader that 
a number of effects make a direct comparison to the RV fre- 
quencies difficult: 1) The RV data are more effected by alias 
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Fig. 9. Schematic of the oscillation spectrum for (3 Gem based on the RV measurements (top) and MOST photometry 
(bottom). 



effects than the photometric data. 2) The spectroscopic and 
photometric data were not contemporaneous and the spec- 
troscopic data in particular were taken over a larger time 
span. 3) The photometric amplitudes detected by MOST 
were small and the time span of the MOST observations 
were relatively short for the detection of such modes. A 
more reliable detection of photometric modes may require 
a more extended MOST observing run. In spite of these 
caveats, it is clear that several of the same oscillation fre- 
quencies are seen in both data sets. 



Figure [S] shows the phase diagrams of the individual 
photometric modes. As in the case for the RV figure the 
contribution of all the other frequencies were first removed 
before phasing the data. For clarity we averaged the data in 
bins of 0.05 in phase. In some cases the frequency and the 
sampling resulted in a natural "clumping" of points which 
dictated a different binning of the data (e.g. / p n). 
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Mode v v Amp FAP 
(cd" 1 ) (/iHz) (PPT) 

/ p i 2.11 ± 0.01 24.41 ± 0.12 0.180 ± 0.013 <10 _l5 

f p2 0.38 ± 0.01 4.40 ± 0.14 0.167 ± 0.012 f»10^ 4 

f p3 6.34 ± 0.01 73.35 ± 0.12 0.142 ± 0.010 f»10~ 5 

f p4 7.50 ± 0.01 86.77 ± 0.12 0.123 ± 0.009 f^lO" 5 

f pB 8.79 ± 0.01 101.70 ± 0.13 0.106 ± 0.011 0.002 

f p6 11.41 ± 0.02 132.01 ± 0.20 0.078 ± 0.017 0.002 

f p7 8.18 ± 0.01 94.64 ± 0.15 0.120 ± 0.013 0.002 

fps 3.39 ± 0.01 39.22 ± 0.16 0.115 ± 0.014 0.002 

f p9 4.09 ± 0.01 47.32 ± 0.11 0.122 ± 0.010 0.002 

/ p io 0.88 ± 0.01 10.18 ± 0.14 0.138 ± 0.012 0.007 

/ p n 2.36 ± 0.03 27.30 ± 0.29 0.132 ± 0.025 0.003 

/pi2 7.04 ± 0.02 81.45 ± 0.23 0.075 ± 0.020 0.002 

/pis 1.69 ± 0.02 19.55 ± 0.26 0.089 ± 0.023 0.027 

/pi4 1.27 ± 0.02 14.69 ± 0.20 0.093 ± 0.017 0.010 

/pis 17.84 ± 0.02 206.40± 0.20 0.072 ± 0.017 0.021 

,/pis 9.57 ± 0.02 110.72 ± 0.23 0.058 ± 0.020 0.010 

Table 4. Pulsation frequencies derived from MOST pho- 
tometry. Amplitudes are in parts per thousand (PPT). 
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Fig. 8. Phase diagrams showing the residual photomet- 
ric measurements phased to the frequencies found in the 
Fourier analysis. For each panel the contribution of all other 
frequencies were subtracted prior to phasing the data. Data 
were phased in phase bins of 0.05 except for f p n where a 
bin of 0.1 was used. 



5. The p-mode Oscillation Spectrum 

The top panel of Figure [9] shows the schematic of the RV 
oscillation spectrum (amplitudes and frequencies) as solid 
lines. The frequency range 70 - 100 /iHz seems to show 
modes that are equally-spaced in frequency with a spacing 
of 10 /iHz However, not all the detected modes seem to 
follow this spacing. We suspect that this may be either due 
to "missing modes" that were not detected either because 
they have too low an amplitude, were simply not present, 
or could not be detected due to the sampling of our data. 



Since we are most likely dealing with p-mode oscillations 
it is reasonable to expect that the modes are evenly spaced 
in frequency. We tried to find a mean frequency spacing 
that could fit the observed spacing. We took the dominant 
mode at 87.24 /iHz as our "origin" and calculated a set of 
predicted frequencies according to v = vq ± nAu, where vq 
= 87.24 /iHz and n an integer. Thus the modes on either 
side of vq will have n — ±1, ±2, .., etc. A value of n was 
chosen so that the predicted mode was as near as possible 
in frequency to the observed modes. A frequency spacing 
Av was then found that minimized the differences between 
the observed and predicted frequencies. 

The dashed lines show the location of the predicted fre- 
quencies using a frequency spacing of Av = 7.14 //Hz. Most 
modes fit extremely well on this frequency spacing with 
the possible exception of the lowest frequency modes. This 
is not surprising since only high order modes should have 
an evenly-spaced frequency separation. The rms scatter of 
the observed minus fitted frequencies (O-C) is 0.38 /iHz. 
Approximately 11 modes were used to derive the frequency 
spacing, so we adopt a value of 0.38/-\/T0 /iHz = 0.12 /iHz 
as the error in the frequency spacing determination. 

The bottom panel of Figure |H] shows the schematic of the 
photometric oscillation spectrum derived from the MOST 
photometry. Again we tried to find a mean frequency sep- 
aration using the same procedure as the RV data. This 
resulted in a frequency spacing for the photometric modes 
with v > 60 /iHz of Av — 7.07 ± 0.32 /iHz consistent with 
the value determined from the RV data. The grid shows 
the location of the predicted frequencies found from the 
RV measurements. 

6. The Mass of /3 Gem 

The large separation of p-modes are given by 

v nit « Av {n + 1/2 + Sv) 

where n and I are the order and degree of the modes, and 
8v is the small spacing that is related to the stellar struc- 
ture (Tassoul 1980). If the modes are non-radial then the 
observed frequency spacing is just one- half of the large spac- 
ing. If the modes are radial (£ = 0) then the large spacing 
is the observed frequency spacing. 

This large separation is related to the mean density of 
the star by 

Av = 135(M/M Q )°- 5 (i?/i? )- 1 - 5 /iHz 

Assuming that these are all radial modes (£ = 0), the spac- 
ing of 7.14 ± 0.12 /iHz and the measured radius of R = 8.8 
± 0.1 R Q results in a stellar mass of M = 1.91 ± 0.09 M Q . 
Using the spacing found in the photometric data results in 
M = 1.88 ± 0.09 Mq. A more precise mass determination 
requires asteroseismic modeling of the observed frequency 
spectrum. 

According to the scaling relationships the frequency of 
maximum power, v max , is also related to the mass and ra- 
dius of the star. From Kjeldsen & Bedding (1995): 

iW = M/R 2 ^/T cS /5777K 3.05 /iHz 

where M, and R are in solar units. So, in principle one can 
determine the mass and radius of the star from v max and 
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Av alone. For the case of /3 Gem, v max = 0.087 mHz (ignor- 
ing low frequency modes) and Av = 7.14 /iHz. This results 
in R = 9.3 Rq and 2.25 M , in good agreement with the 
interferometric radius and subsequent mass determination. 

We should note that these scaling relations were derived 
from solar type stars, but seem to hold also for evolved 
red giant stars (Kallinger et al. 2009; Stello et al. 2009). 
Also, recent results from the Kepler Mission indicate that 
K giants can have mixed 1 = 1 modes (Mosser et al. 2012) 
which can shift the observed frequencies. We cannot say 
from our data whether (3 Gem has such modes or whether 
we could have detected them with our RV measurements. If 
present they may affect the accuracy of frequency spacing 
determination. 

7. Discussion 

Our long time series of precise RV measurements for /3 Gem 
reveals the clear presence of p-mode oscillations in this star. 
These frequencies fall on a grid of equally-spaced frequency 
on an interval of Av = 7.14 ± 0.2 /iHz. Photometric mea- 
surements taken with the MOST micro-satellite also show 
a comparable frequency spacing. Most of the photomet- 
ric mode match up well with those detected in the spec- 
troscopy. 

It seems that these modes must be high overtone radial 
modes (£ = 0). Nonradial modes would have a large spacing 
of twice the observed frequency spacing or ss 14 /iHz. This 
implies a stellar mass M « 8 Mq which would be greatly 
at odds with stellar evolutionary tracks. 

There is evidence that K giant stars can exhibit both 
radial and nonradial modes. Hatzes & Cochran (1994) sug- 
gested that the short-period RV variations in a Boo were 
consistent with fundamental or low-overtone radial modes. 
Buzasi et al. (2000) reported the detection of 10 modes 
in a UMa using the guide camera of the failed WIRE 
satellite. These were spaced on frequency interval of 2.94 
/iHz which was consistent with overtone radial pulsations. 
Radial modes have also been claimed for e Oph (De Ridder 
et al. 2006) using RV measurements. Barban al. (2007) ob- 
served this star photometrically with MOST and detected 
7 equidistant peaks that were interpreted as radial modes. 
However, this result was not without dispute. Kallinger et 
al. (2008) re-examined the available data on e Oph and 
argued for the presence of 21 modes. They were able to 
match these modes with a model that included both radial 
and non-radial modes. The unequivocal detection of both 
radial and non-radial modes finally came from light curves 
from the CoRoT space telescope. De Ridder et al. (2009) 
examined light curves from giant stars from the first long 
run of 150 days (LRcOl) from CoRoT and found strong ev- 
idence for the presence of radial and nonradial modes, par- 
ticularly in the star CoRoT-101034881. However, Baudin 
et al. (2012) used 60-d of CoRoT observations of the giant 
star HD50890 and found oscillations between 10-20 /iHz 
that were consistent with a regular spacing of 1.7 /iHz con- 
sistent with only radial modes. So it may be that some 
giant stars can oscillate in predominantly radial modes and 
others with both modes. (Of course the stars with pure ra- 
dial pulsations may also possess non-radial modes that are 
simply not detected.) 

There is some evidence for nonradial modes, at least in 
the photometric data. For instance, f p u, f p i5 and possibly 
/ P 9 may lie on the half interval of the large spacing derived 



from what we think are radial modes. Nonradial modes may 
well exist in (3 Gem, but possibly were not detected because 
of their low RV amplitudes and/or because of our limited 
time coverage. Furthermore, until the modes are actually 
identified some of the ones shown in Figure [5] may in fact 
be nonradial modes. 

Evidence for the type of modes (non-radial versus ra- 
dial) can also come from the ratio of radial velocity to 
photometric amplitudes (2K/AV). Radial pulsators such 
as Cepheids tend to have 2K/AV « 50-100 kms -1 mag -1 
while nonradial pulsators tend to have smaller values. Table 
5 lists the frequencies in the RV data for which we think 
we have found a counterpart in the photometric data. The 
2K/AV ratio is listed in the last column. All but one mode 
have 2K/AV in the range of 57-103 consistent with radial 
pulsators. One mode (fvio/fpu) shows a 2K/AV ratio con- 
sistent with nonradial modes. We caution the reader that 
because the photometric and spectroscopic data were not 
simultaneous we cannot be sure that the RV amplitude for 
a mode was the same during the time of the photomet- 
ric measurements. A more accurate determination of the 
radial velocity to photometric amplitude would require si- 
multaneous measurements. However, at face value both the 
frequency spacing of the modes and the values 2K/AV ar- 
gue in favor of radial modes. 

It is not entirely clear what the nature of the low fre- 
quency modes seen in both the spectroscopic and photo- 
metric data is, but low order radial pulsations may be a 
prime candidate. The expected frequencies of the funda- 
mental (F) and first 2 harmonics (1H, 2H) for radial pul- 
sations are 1.5, 2 and 2.5 cd" 1 (P = 0.67, 0.5, and 0.4 d) 
using the pulsation constants of Cox, King, & Stellingwerf 
(1972). These are close to f v3 (2H) and f vl4 (F) in the RV 
data and f pl (2H) and f pU (F) in the MOST data. Using 
the pulsation constants from Xiong & Deng (2007) which 
may be more appropriate for cool giants results in a funda- 
mental and first harmonic frequency of 1.67 and 3.22 c d _1 , 
respectively. These also close to the frequencies of /„i4 and 
f V 9 in the RV data and / p i3 and f p % in the photometry. If 
these are indeed low order radial modes then it seems that 
f3 Gem can excite radial modes over a wide range of orders. 
We note that the low frequency modes f v \ and f v § in the 
RV and f P 2 and f p -j in the photometry cannot readily be 
attributed to radial pulsations. However, a detailed pulsa- 
tional analysis using appropriate stellar models for f3 Gem 
arc needed to identify properly all detected modes. 



VRV 


Amp/iv- 


VPhot 


Ampph ot 


2K/Am 


(/iHz) 


(ms- 1 ) 


(/iHz) 


(milli-mag) 


(kms - mag -1 ) 


80.49 


2.00 


81.45 


0.075 


53 ± 14 


87.24 


4.31 


86.77 


0.123 


70 ± 5 


93.02 


4.00 


94.64 


0.112 


71 ± 6 


10.03 


5.39 


10.18 


0.138 


78 ± 8 


71.58 


1.05 


73.35 


0.140 


15 ± 7 


130.88 


1.50 


132.01 


0.078 


38 ± 14 



Table 5. Ratio of RV to photometric amplitude of modes 
found in j3 Gem. 



Besides radial versus nonradial modes, the lifetime of 
the modes in K giant stars has also been the subject of 
discussion in the literature. By fitting Lorentz profiles to 
the power spectrum of e Oph Barban et al. (2007) derived 
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mode lifetimes of 2.7 days. However, Kallinger et al. (2008) 
suggested lifetimes of 10 to 20 days and argued that many 
of the peaks in the broad Lorentz profiles may have actually 
been independent modes. Photometric data on giant stars 
obtained by CoRoT indeed confirmed that some K giant 
stars can have mode lifetimes of tens of days (De Riddcr 
et al. 2009), if not longer. Our excellent fitting of our RV 
data for /3 Gem that spans 120 days also suggest the mode 
lifetimes, at least for this star, are relatively long, or are 
re-excited on a regular basis. 

The main goal of our analysis was to derive the p-mode 
frequency spacing of /3 Gem and to combine this with the 
interfcrometric stellar radius in order to determine a model 
independent stellar mass. Our resulting stellar mass of 1.91 
± 0.09 Mq is in excellent agreement with the value of 1.96 
± 0.19 M determined from evolutionary tracks. However, 
we should note that the mass determination from evolution- 
ary tracks depends critically on the derived stellar param- 
eters. For example, if we use the T e g and [Fe/H] of Allcnde 
Prieto et al. (2004) we obtain a stellar mass of M = 1.47 ± 
0.45 Mq. Thus the actual error in the mass determination 
from evolutionary tracks is probably much larger than the 
formal errors using a given set of stellar parameters. The 
stellar mass derived from the frequency spacing is thus a 
better determined value since it is independent of T c s and 
abundance. /? Gem is thus one of the few planet hosting 
stars for which we have an accurate measurement of the 
stellar mass. The combination of stellar oscillations and in- 
terferometric radius determinations can provide a powerful 
tool for verifying stellar evolution codes using nearby, bright 
giant stars and to obtain more accurate stellar masses for 
planet hosting giant stars. 

The questions regarding whether /3 Gem oscillates only 
in radial modes, whether it can also excite nonradial modes, 
or what the actual mode lifetimes are all remain open. 
Unfortunately, our RV data set has several long data gaps 
and the MOST data length is too short. Long, uninter- 
rupted observations are required for this and this can only 
be done photometrically from space, or spectroscopically 
from multi-site campaigns or with a dedicated network 
of telescopes. This star is an ideal target for spectro- 
scopic studies with RV measurements with the proposed 
SONG network of telescopes (Grundahl et al. 2009) or with 
high precision photometry using the micro-satellite BRITE- 
Constellation (Weiss et al. 2008). 
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